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, The generalized Henon-Heiles system has been considered. In two nonintegrable cases with 

the help of the Painleve test new special solutions have been found as Laurent series, de- 
pending on three parameters. The obtained series converge in some ring. One of parameters 
determines the singularity point location, other parameters determine coefficients of series. 
For some values of these parameters the obtained Laurent series coincide with the Laurent 
■ series of the known exact solutions. The Painleve test can be used not only to construct 

OO ' local solutions as the Laurent series but also to find elliptic solutions. 

o 

(N 

g : 1 THE PAINLEVE PROPERTY AND INTEGRABILITY 



A Hamiltonian system in a 2s-dimensional phase space is called completely integrahle (Liouville 
integrable) if it possesses s independent integrals which commute with respect to the associated 
Poisson bracket. When this is the case, the equations of motion are separable (at least, in 
principal) and solutions can be obtained by the method of quadratures. 

When some mechanical problem is studied, time is assumed to be real, whereas the integra- 
bility of motion equations depends on the behavior of their solutions as functions of complex 
time. S.V. Kovalevskaya was the first, who proposed ^ to interpret time as a complex vari- 
\^ ' able and to require that the mechanical-problem solutions have to be single-valued functions 



X 



meromorphic in the entire complex plane. This idea led S.V. Kovalevskaya to a remarkable 
result: a new integrable case (nowadays known as the Kovalevskaya's case) for the motion of 
a heavy rigid body about a fixed point was discovered P (see also jHEl)- The Kovalevskaya's 
result demonstrated that the analytic theory of differential equations can be fruitfully applied 
to mechanical and physical problems. The important stage of development of this theory was 
the Painleve classification of ordinary differential equations (ODE's) with respect to the types 
of singularities of their solutions j^. 

Let us formulate the Painleve property for ODE's. Solutions of a system of ODE's are 
regarded as analytic functions, may be with isolated singularity points |S1 E]. A singularity 
point of a solution is said critical (as opposed to noncritical) if the solution is multivalued 
(single-valued) in its neighborhood and movable if its location depends on initial conditions. 
The general solution of an ODE of order N is the set of all solutions mentioned in the existence 
theorem of Cauchy, i.e. determined by the initial values. It depends on arbitrary independent 
constants. A special solution is any solution obtained from the general solution by giving values 
to the arbitrary constants. A singular solution is any solution which is not special, i.e. which 
does not belong to the general solution. 



Definition 1. A system of ODE's has the Painleve property if its general solution has no 
movable critical singularity point [1]. 
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Investigations of many dynamical systems ||2| show that systems with the Painleve prop- 
erty are completely integrable. Arguments, which clarify the connection between the Painleve 
analysis and the existence of motion integrals, are presented in 131^]. At the same time the 
integrability of an arbitrary system with the Painleve property has yet to be proved. There is 
not an algorithm for construction of the additional integral by the Painleve analysis. It is easy 
to give an example of an integrable system without the Painleve property jTH]. The system with 
the Hamiltonian H = + f{x), where f{x) is a polynomial which power is not lower than 
five, is trivially integrable, but its general solution is not a meromorphic function. 

The study of complex-time singularities is a useful tool for the analysis of not only integrable 
systems, but also chaotic dynamics The Painleve analysis can be connected with the normal 
form theory 1 1 2J . 

The Painleve test is any algorithm, which checks some necessary conditions for a differ- 
ential equation to have the Painleve property. The original algorithm, developed by P. Painleve 
and used by him to find all the second order ODE's with the Painleve property, is known as the 
a-method. The method of S.V. Kovalevskaya is not as general as the a-method, but much 
more simple. The remarkable property of this test is that it can be checked in a finite number of 
steps. This test can only detect the occurrence of logarithmic and algebraic branch points. Up 
to the present there is no general finite algorithmic method to detect the occurrence of essential 
singularities. Different variants of the Painleve test are compared in jl3| R. Conte paper]. 

Developing the Kovalevskaya method .1 further, M.J. Ablowitz, A. Ramani and H. Segur 
constructed a new algorithm of the Painleve test for ODE's |T1]. They also were the first to 
point out the connection between the nonlinear partial differential equations (PDE's), which 
are solvable by the inverse scattering transform method, and ODE's with the Painleve property. 
Subsequently the Painleve property for PDE's was defined and the corresponding Painleve test 
(the WTC procedure) was constructed [KJ (see also ^1 EI ) • With the help of this test it 
has been found, that all PDE's, which are solvable by the inverse scattering transforms, have 
the Painleve property, may be, after some change of variables. For many integrable PDE's, 
for example, the Korteweg-de-Vries equation the Backlund transformations and the Lax 
representations result from the WTC procedure jl6lll8j . For certain nonintegrable PDE's special 
solutions were constructed using this algorithm I 19| 120] . 

The algorithm for finding special solutions for ODE's in the form of finite expansions with 
respect to powers of an unknown function has been constructed in [21\ I22j. This function and 
coefficients have to satisfy some system of ODE's, often more simple than an initial one. This 
method has been used to obtain exact special solutions for some nonintegrable ODE's |23| I24j . 
With the help of the perturbative Painleve test ^21 a four-parameter generalization of an exact 
three-parameter solution of the Bianchi IX cosmological model has been constructed j25) . 



In the 1960s the models of the star motion in an axial-symmetric and time-independent potentials 
have been developed to show either existence or absence of the third integral for some polynomial 
potentials. Due to the symmetry of the potential the considered system is equivalent to two- 
dimensional one. To clarify the question about the existence of the third integral Henon and 
Heiles [2n] considered the behavior of numerically integrated trajectories. Emphasizing that 
their choice does not proceed from experimental data, they have proposed the Hamiltonian 



because on the one hand, it is analytically simple; this makes the numerical computations of 
trajectories easy; on the other hand, it is sufficiently complicated to give trajectories which are 



2 THE HENON-HEILES HAMILTONIAN 
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far from trivial. Indeed, for low energies the Henon-Heiles system appears to be integrable, 
in so much as trajectories (numerically integrated) always lay on well-defined two-dimensional 
surfaces. On the other hand, for high energies many of these integral surfaces are destroyed, it 
points on absence of the third integral. 

The generalized Henon-Heiles system is described by the Hamiltonian: 

H = ^(^xf + y^ + Xx^ + y^^+x^y-jy^ (1) 

and the corresponding system of the motion equations: 

f xtt = -Ax - 2xy, 

2 2 

I ytt = -y- X +Cy , 

where xu = ^ and yu = A and C are numerical parameters. 

The generalized Henon-Heiles system is a model not only actively investigated by various 
mathematical methods (see and references therein), but also widely used in astronomy and 
physics, in particular, in gravitation [281 129j . The models, described by the Hamiltonian 
with some additional nonpolynomial terms, are actively studied ^Hl 1^ OH as well. 

Due to the Painleve analysis the following integrable cases of have been found: 

(i) C = -l, A=l, 

(ii) C = —6, A is an arbitrary number, 

(iii) C = -16, A = 1/16. 

The general solutions in the analytic form are known only in the integrable cases j31l I32j , in 
other cases not only four-, but even three-parameter exact solutions have yet to be found. In 
nonintegrable cases local four-parameter solutions as converging psi-series solutions have been 
found 33j for all values of the parameter C, except C = —2. The Ablowitz-Ramani-Segur 
algorithm of the Painleve test appears very useful to find such values of parameter at which 
three-parameter solutions can be expanded in formal Laurent series and to construct these local 
solutions. The knowledge of local solutions can be used to find solutions in the analytical form. 

Let us assume that the behavior of solutions in a sufficiently small neighborhood of the 
singularity point is algebraic and x and y tend to infinity as 

X = aa{t — to)"' and y = bjs^t — to)^ , 

where a, /3, and 6^ are some constants. 

The Painleve test gives all information about behavior of solutions in the neighborhood of the 
singularity point (see, for example, [7])- There exist two possible variants of dominant behavior 
and resonance structure of solutions of the generalized Henon-Heiles system (3 ES] : 



Case 1: 


Case 2: (/3 < 3?e(a)) 


a = -2, 
/3 = -2, 

a„ = ±3V2 + C, 

bp = -3, 

1 5 Vl-24{1+C) 5 1 Vl-24{1+C) 
/ — 1, U, 2 2 ' 2 ' 2 


_ l±Vl-48/C 

/3 = -2, 

fl« = (an arbitrary number), 
h - ^ 

r = -l, 0, 6, T^l-f 



The values of r denote resonances: r = — 1 corresponds to arbitrary parameter ioj = (in 
the Case 2) corresponds to arbitrary parameter c^. Other values of r determine powers of t, to 
be exact, for x and t^^^ for y, at which new arbitrary parameters can appear. 
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For integrability of system ((2j all values of r have to be integer and all systems with zero 
determinants have to have solutions at any values of included in them free parameters. It is 
possible only in the integrable cases (i) — (iii). 

Those values of C, at which r are integer numbers either only in the Case 1 or only in the Case 
2, are of interest for search of special three-parameter solutions. Those cases where an additional 
negative resonance is present likely correspond to singular, rather the general, solutions 0. 

Let's consider all cases, when there exist special (no singular) solutions, representable as a 
three-parameter Laurent series (may be, multiplied on y/t — to). From the requirement that all 
values of r but one are integer and nonnegative numbers we obtain the following values of C: 

C = -1,C = -4/3 (the Case 1), C = -16/5, C = -6, C = -16 (the Case 2, a = IzVlE^^ 
and C = —2, when two types of singular behaviour coincide. 

At C = —2 (in the Case 1) aa = 0. This is the consequence of the fact that, contrary to 
our assumption, the behaviour of the solution in the neighborhood of a singularity point is not 
algebraic, because its dominant term includes logarithm [Jj. At C = —6 and any value of A the 
exact four-parameter solutions are known. In cases C = —1 and C = — 16 the substitution of 
unknown function as Laurent series gives the equations in A: accordingly A = 1 and A = 1/16, 
hence, in nonintegrable cases special three-parameter local solutions have to include logarithmic 
terms. Single-valued three-parameter special solutions can exist only in two nonintegrable cases, 
at C = -16/5 and at C = -4/3. 

3 NEW SOLUTIONS 

Let us consider the Henon-Heiles system with C = —16/5. In the Case 1 some values of r are 
not rational. To find special three-parameter solutions we consider the Case 2. In this case 
a = —3/2 and r = —1, 0, 4, 6, hence, in the neighborhood of the singularity point we have 
to seek x in such a form that can be expanded into Laurent series, beginning with (t — to)~^- 
Let = 0, substituting 

(oo \ oo 

cit"2 -h ^ ajt^ and y = - — i"^ + ^ hjt^ (3) 

i=-i / i=-i 

in 0, we obtain the following sequence of linear system in and h^: 

fc-i 



i=-i 

(4) 

(^{k - l)k - I2^bk = - bk-2 - ^ ajak-j-3 " y bjbk-j-2- 

i=-2 j=-i 

The determinants of the systems (01) corresponding to A; = 2 and k = A are equal to zero. To 
determine 02 and 62 we have the following system: 

ci (^557056c? + (15552000A - 4860000) cf + 864000000^2 + 

+ lOSOOOOOOA^ - 67500000A + 10546875^ = 0, (5) 
818176c? + (^15660000A - 489375o)cf - 8IOOOOOOO62 - 6328125 = 0. 

It is easy to see that this system contains no terms proportional to 02, therefore, 02 is the 
new constant of integration. We discard the solution with = and obtain the system in c\ 
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and 62- System (jSJ has solutions only if 

^ 1125 (^525 - 1680A ± 4y/35(2048A2 - 1280A + 387)) 
~ 167552 ■ 

We obtain new constant of integration 02, but we must fix c^, so number of constants of 
integration is equal to 2. It is easy to verify that 64 is an arbitrary parameter, because the 
corresponding system is equivalent to one linear equation. System (j2| is invariant under exchange 
X to —X, so, we obtain four different local solutions which depend on three parameters, namely 
to, 02 and 64. With the help of some computer algebra system, for example, REDUCE |34j . 
these solutions can be obtained with arbitrary accuracy. For the case A = 1/9 the obtained 
Laurent series are presented in |35j . 

At C = —4/3 the situation is similar. In the Case 1 we have r = —1, 1, 4, 6. Substituting 
00 00 
X = Vet^^ + dkt^ and y = -St'^ + ^ fkt^ (6) 
fc=-i fc=-i 

in system we receive a sequence of linear systems in dk and fk- 

k-l 

\k - l)k - 6)4 + 2\/6A = -Xdk-2 - 2 <i,/fc-j-2, 



fc-i ^ fc-i 



(7) 



2V6dk + Uk - l)k - 8j/fe = - fk-2 - Yj djdk-j-2 - fjfk-j-2- 



3 

i=~i i=-i 



The determinants of the systems corresponding to A; = —1, 2, 4 are equal to zero. The first 
system (A; = —1) always has infinite number of solutions and /_i is a parameter. We have to fix 
this parameter to solve the system corresponding to A; = 2. This system has solutions only if 



,2 105 - 140A ± v^7(1216A2 - 1824A + 783) 

/ 1 = or /_i = 0. 

•'-^ 385 ■' 

At k = A system Q is reduced to one equation. Thus, at C = —4/3 we have five three- 
parameter (to, /2 and fi) solutions. 

The convergence of all the Laurent series solutions on some real time interval have been proved 
in j33j . For the obtained solutions it is easy to find conditions, at which the series converge at 
< |t| ^ 1 — e, where e is any positive number. Our series converge in the above-mentioned 
ring, if 3N and 3M such that Vn > iV |a„| ^ M and ^ M. Let |a„| ^ M and |6„| ^ M for 
all —1 < n < k, then (in the case C = —16/5) from (jlj) we obtain: 

, . 2M(A; + l) + |A|+2|cil ^ . 21MA; + 26M + 5 

^ \^^\ ^' ^ 5\k^-k-12\ 

It is easy to see that there exists such that if |a„| ^ M and ^ M for —1 ^ n ^ N, 
then |a„| ^ M and ^ M for — 1 ^ n < cxd. So one can prove the convergence, analyzing 
values of a finite number of the first coefficients of series. For C = —4/3 it is easy to obtain the 
analogous result. 



4 Global single- valued solutions 

We have found some local three-parameter solutions. To seek the global single- valued solutions 
we transform system Q into the equivalent fourth order equation (351 13 7j : 

20c 

ytttt = {2C - ^)ytty - (4A + l)ytt + 2(C + l)yl + + (4C7A - Q)y^ - 4Ay - 4F, (8) 
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where H is the energy of the system. We note, that the value of H depends on initial data, that is 
to say that from five parameters: y(io)) yt{to), Vttito)-, Utttito) and H only four are independent. 

There are some reasons to seek three-parameter solutions of eq. © in terms of elliptic func- 
tions. In 1999 E.I. Timoshkova |SZI found that the general solution of the following equation: 



with some values of constants A, B, C, V, Q and <S, is one-parameter solution of eq. © in each 
of two above-mentioned nonintegrable cases (C = —4/3 or C = —16/5, A is an arbitrary 
number). If ^ = and if = we obtain the well-known solutions in terms of the Weierstrass 
elliptic function. Solutions with Q ^ {) or £ are derived only at P = 0, therefore, substitution 
y{t) = Q{tf gives: 



Two-parameter solutions y{t) = g^t)"^ + -Po; where Pq is an arbitrary constant and Q{t) satisfy 
equation ©, have been obtained in [HSllSni- These solutions are the following elliptic functions 



where p(t — Iq) is the Weierstrass elliptic function, a, 6, c and d are some constants. The 
parameter Pq defines the energy of the system. There exist two different elliptic solutions for 
each possible pair of values of C and A. 

Let us consider the three-parameter solutions which were obtained at C = —4/3. If we 
choose /_i = 0, then we obtain the solution which generalizes the known two-parameter solution 
in terms of Weierstrass elliptic functions. Other solutions generalize two-parameter solutions, 
obtained in [2^1 • The coefficient /_i is a residue of y. The sum of residues of elliptic function 
in its parallelogram of periods has to be zero , hence, two local solutions with opposite signs 
of /_i correspond to one global elliptic solution. The obtained local three-parameter solutions 
generalize the Laurent series of the two-parameter elliptic solutions in the form (jlUj) . 

For C = —16/5 we obtained four local solutions, which generalize two global elliptic solutions 
in the form (|1U() . So, each obtained local three-parameter solution generalize the Laurent series of 
some two-parameter elliptic solution and we can assume that an unknown global there-parameter 
solutions are elliptic functions. 

Of course, solutions, which are single-valued in the neighborhood of one singularity point, 
can be multivalued in the neighborhood of another singularity point. So, we can only assume 
that global three-parameter solutions are single- valued. If we assume this and moreover that 
these solutions are elliptic functions (or some degenerations of them), then we can seek them as 
solutions of some polynomial first order equations. The classical theorem, which was established 
by Briot and Bouquet j41j . proves that if the general solution of the autonomous polynomial 
ODE is single- valued, then this solution is either an elliptic function, or a rational function of 
e^^ , 7 being some constant, or a rational function of x. Note that the third case is a degeneracy 
of the second one, which in its turn is a degeneracy of the first one. At the same time, there 
exist elementary functions, for example, the function f[t) = t + sin(t), which are not solutions 
of any first order polynomial ODE. 

A second result, of immediate practical use, is due to P. Painleve 4^. He has proved that 
if the general solution of an autonomous polynomial ODE is single- valued, then the necessary 
form of this ODE is 



2/2 = ^ ^y2 ^Cy + V + Gy^l^ + ,Sy3/2 




(9) 




[ad — 6c = 1) 



(10) 



m 2m— 2k 




(11) 



A:=0 j=0 
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in which m is a positive integer number and hjk are constants. 

In 2003 R. Conte and M. Musette have proposed a new method to find elliptic solutions |42j . 
This method is based on the Painleve test and uses the Laurent series expansion to find the 
analytic form of elliptic solutions. Rathen than to substitute the first order equation ^11^ into 
equation (jSJ one can substitute the found Laurent series solutions of equation Q, for example, 
either solution or solution ®, into equation (|1H) and obtain a linear system in hjk- This 
method is more powerful than the traditional method and allows in principal to find all elliptic 
solutions. I hope that the use of this method allows to find the three-parameter elliptic solutions. 



5 CONCLUSION 

Using the Painleve analysis we have found local special solutions in two nonintegrable cases of 
the generalized Henon-Heiles system (C = —16/5 and C = —4/3). These solutions are the 
converging Laurent series, depending on three parameters. For some values of these parameters 
the obtained solutions coincide with the known exact periodic solutions. There are no obstacles 
to exist three-parameter single-valued solutions, so, the probability of finding of exact three- 
parameter solutions, which generalize the solutions obtained in [371 139j . is high. 

The author is grateful to F. Calogero, R. Conte, V. F. Edneral, A. K. Pogrebkov and 
E. L Timoshkova for valuable discussions. This work has been supported by the Russian Foun- 
dation for Basic Research under grants 00-15-96560, 00-15-96577 and NSh-1685.2003.2 and by 
the grant of the scientific Program "Universities of Russia" UR. 02. 03. 002. 
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